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Q" 1 Abstract 

^Jj This paper deals with the dynamics of scalar field thin shell in the 

Reissner-Nordstrom geometry. The Israel junction conditions between 
Reissner-Nordstrom spacetimes are derived, which lead to the equa- 
tion of motion of scalar field shell and Klien-Gordon equation. These 
equations are solved numerically by taking scalar field model with the 
quadratic scalar potential. It is found that solution represents the ex- 
panding and collapsing scalar field shell. For the better understanding 
of this problem, we investigate the case of massless scalar field (by tak- 
ing the scalar field potential zero). Also, we evaluate the scalar field 
potential when p is an explicit function of R. We conclude that both 
massless as well as massive scalar field shell can expand to infinity at 
constant rate or collapse to zero size forming a curvature singularity 
or bounce under suitable conditions. 

Keywords: Gravitational collapse, Israel junction conditions; Scalar field. 
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1 Introduction 

The idea of the geon (electromagnetic-gravitational entity) by Wheeler et al. 
[HE] leads to investigation of scalar field in general relativity (GR). In GR, 
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the scalar field appears in the low energy limit of string theory [3J. Still, 
there is no observational evidence about the existence of such particles that 
are associated with the scalar field. However, the study of gravitational col- 
lapse of compact objects in the scalar-tensor theories imply that scalar field 
might be the source of scalar gravitational waves that can be detected by 
the advanced detectors Further, the observational facts from the binary 
pulsar may include or exclude the presence of scalar field in GR [3] . In astro- 
physical context, boson stars are such compact objects that are composed of 
scalar field. Probably, such stars were created in the early universe as Higgs 
particles condensate [6] . Recently, it has been proposed that boson stars are 
the candidates for dark matter [7J. Also, some remarkable similarities are 
suggested between neutron and boson stars. 

Dynamics of scalar fields has been the subject of particular interest in 
both cosmological as well as astrophysical situations. In cosmological scenar- 
ios, scalar fields have a great attraction because such fields act as an effective 
cosmological constant in deriving the inflation. The nature of spacetime 
singularity for massless scalar field has been investigated in recent years on 
spherical collapse models [SJ-[I5]. Bhattacharya et al. [IS] discussed the grav- 
itational collapse of a minimally coupled massless scalar field and examined 
the possibility of existence of nonsingular models where collapse could be 
freezed under suitable conditions. All these studies investigate the dynamics 
of massless scalar field using the exact solutions of the field equations. Kaup 
[5] was among the pioneers to study complex massive scalar field configura- 
tion. Ruffini and Bonazzola [TTJ explored spherically symmetric system and 
determined the equilibrium conditions for boson stars solutions. 

There is another formalism to study the dynamics of the matter field re- 
ferred to as " thin shell formalism" developed by Israel [IB] . This is one of 
the exactly solvable formulation in GR which is widely used to understand 
gravitational collapse and other cosmological as well as astrophysical pro- 
cesses. It involves a discontinuity of the extrinsic curvatures of the interior 
and exterior regions across a boundary surface. The jump in the extrinsic 
curvature across the boundary surface is caused by a thin layer of matter. 
In this formulation, the set of equations leads to equations of motion (much 
simpler) corresponding to the field equations, involving curvature on one side 
and matter on the other side. The solution of these equations provides full 
understanding of the dynamics of system. In the relativistic astrophysics, 
the thin shell equations help to study the properties of the compact objects. 

Pereira and Wang [19] studied gravitational collapse of the cylindrical 
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shell made of counter rotating dust particles by using the Israel thin shell 
formalism. Sharif et al. [20j-[22j have investigated plane and spherically 
symmetric gravitational collapse by using this formulation. This approach 
was generalized to thin charged shell without pressure by De La Cruz and 
Israel [23]. Kuchar (23] and Chase [25] treated the charged thin shell prob- 
lem with pressure by using polytropic equation of state. There are a number 
of papers devoted to handle the charged thin shell problems. Boulware [2B] 
studied the time evolution of the charged thin shell and showed that their 
collapse could form a naked singularity if and only if density is negative. 
Farrugia and Hajicek [27] investigated third law of black hole mechanics in 
the Reissner Nordstrom (RN) geometry. Nunez [28J studied the oscillating 
perfect fluid shell. Also, Nunez et al. [29J explored the stability and dynam- 
ical behavior of the real scalar field for the Schwarzschild geometry in single 
null coordinate. 

The purpose of the present paper is to study the application of Israel thin 
shell formalism for the scalar fields in the context of GR. In particular, we 
study the dynamical behavior of scalar field thin shell in charged background. 
To this end, we take interior and exterior regions as RN solutions, the Israel 
thin shell formalism is used to derive equations of motion of the shell for the 
perfect fluid. We then specify that the perfect fluid is generated by scalar 
field which leads to scalar field equations of motion. The solution of these 
equations provide expanding, collapsing or oscillating scalar field shell. The 
plan of the paper is as follows: Equations of motion for scalar field model 
with quadratic potential and their physical interpretation are presented in 
section 2. Sections 3 and 4 deal with dynamics of the massless and massive 
scalar field with arbitrary scalar potential, respectively. We summarize our 
results in the last section. 

2 Dynamical Equations 

In this section, we use Israel thin shell formulation to derive equations of 
motion of the scalar field shell. For this purpose, we take a 3D spacelike 
boundary surface E, which splits the two 4D spherically symmetric space- 
times V + and V~ . The interior and exterior regions l/~and V + , respectively 
are described by the RN metrics given by 

(d s f ± = v± dT 2 - - dR 2 -R 2 (d6 2 + sin 2 6d<J) 2 ), (1) 
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where r]±(R) = 1 — + -^t, M± and Q± are the mass and charge, respec- 
tively. The subscripts + and — represent quantities in exterior and interior 
regions to £, respectively. Further, it is assumed that charge in both regions 
is the same, i.e., Q_ = Q + = Q. The strength of the electric field on the 
shell can be described by the Maxwell field tensor, Ftr = = —F RT . 
The energy-momentum tensor of the electromagnetic field is 

T u{em) = l { _ F »X Fsx + I^F^F^). (2) 
47T 4 

By employing the intrinsic coordinates (r, 9, 0) on £ at R = R(r), the metrics 
(P]) on £ become 

= [v(R) - ^y(^) 2 ] rfT2 " R 2 (r)(d6* + sin 2 9dcj ) 2 ). (3) 

Here T is also a function of r and it is assumed that goo > 0, so that T is a 
timelike coordinate. Also, the induced metric on the boundary surface £ is 
given by 

(ds)\ = dr 2 -a 2 (r)(d9 2 + sin 2 9d(P 2 ). (4) 
The continuity of the first fundamental form gives 

1 d R 

W^)-^)(-^) a ] 5dr = (*■)=' *(T) = a(T) E . (5) 
Now the unit normal to £ in coordinates can be evaluated as 



n /i ± = (- J R(r),T,0,0), (6) 

where dot represents differentiation with respect to r. 

The non-vanishing components of the extrinsic curvature are 

Kt = -^]R? + r ]± , K± = -Ry/fr + ri ± , K^ = K± sin 2 9. (7) 
The surface energy-momentum tensor is defined by 

fy = i{[#y]-7«[Jr]}, ( 8 ) 

ri 

where k is a coupling constant, 7^ is induced metric on £ and 

[K tJ } = Kt J -Kr p [K\=<f[K ij ]. (9) 
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The surface energy-momentum tensor for a perfect fluid is 

<% = {p + p)uiUj ~ Plij, (10) 
where itj = 5°. Using Eqs.([5]), ([H]) and ( TTUj) . we can find 

Inserting the non-zero components of the extrinsic curvature components, 
we get 

(C+ - C-) + \pR = o, (12) 

^(C+ - C-) + ^(C+ - C-) - np = 0, (13) 



where C± = y R 2 + V±- Making use of Eq.(fT2l in ( [TBI , it follows that 

m + pi = 0, (14) 

where m(= pA) and -A(= A7tR 2 (t)) stand for the integrated total energy 
density at some time and area of the shell, respectively. The conservation 
of surface energy- momentum tensor leads to the same equation as Eq. ( 114ft 
and hence this equation is known as energy conservation law on the shell. It 
is mentioned here that this equation can be solved by using the equation of 
state p = kp, and its solution is 

P = Po(f) 2(fc+1) , (15) 

where Rq is the position of the shell at t — to and po is the density of matter 
on the shell at position Rq. Using this value of p in the definition of m, we 
obtain 

R (2k+2) 

From Eq.( ll2p . equation of motion of the shell is given by 

R 2 + V eff (R) = 0, (17) 
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where the effective potential V e ff(R) is 



Notice that we have used k = 8% to derive this equation. 

To see the effects of charge parameter Q on the dynamics of shell, we 
re- write Eq. lfTT]) by using the above equation as follows 

Here +(— ) correspond to expansion (collapse) of shell and m is the same as 
defined after Eq. fll4l) . The term ^ (Coulomb repulsive force) in R (velocity 
of shell with respect to stationary observer) indicates that charge reduces 
velocity of the shell with respect to stationary observer. This velocity also 
depends on position of the observer, whether the observer is located inside or 
outside the shell. Further, Eq.f llSp implies that charge parameter increases 
the effective potential V e ff. In coming sections, we shall see that throughout 
the dynamics of shell (composed of either massless or massive scalar field), 
charge parameter reduces and increases velocity of the shell with respect to 
stationary observer and V e ff, respectively. Thus initially the velocity of the 
shell with respect to stationary observer in the RN background is slower as 
compared to the Schwarzschild case (as shown in Figure 1). We conclude 
that electrostatic repulsive force in RN background tries to balance with the 
gravitational force due to the shell and hence the shell velocity with respect 
to stationary observer is slow in this case as compared to uncharged case. 

In order to study the dynamics of scalar field shell, we specify that the 
perfect fluid is generated by a scalar field. The energy density and the pres- 
sure of the scalar field can be written as 

p = ^[<P,^ u + 2V(<j ) )], p=±[<l>, u <l>> v -2V(<l>)], (20) 

where V((f>) = m 2 (p 2 is the potential term which contributes to provide mass 
of the scalar field. We note that the scalar field will be massless in the absence 
of such term. Using Eq. fpZOl) . we can write the energy-momentum tensor of 
the scalar field as follows 

= V^V> - 7o - [^(V0) 2 - V(<j>)\. (21) 
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Figure 1: Behavior of the shell velocity with respect to stationary observer, 
when M + = 1, M_ = 0, k = po = Ro = 1 and Q — 1. The upper and 
lower curves correspond to uncharged and charged cases, respectively. It is 
clear that initially velocity in the charged case is less than the uncharged 
case. Velocities in both cases match for larger values R, as term ^ becomes 
negligible for larger values R 




Figure 2: The left graph is the behavior of the shell radius R while the right 
graph is the behavior of scalar field. Both these graphs have been plotted by 
using M+ = 1, m = 1, M_ = 0, Q = 1, 0(1) = 0.19 and 0(1) = R(l) = 1. 



Since the induced metric (4) depends only on r, so also depends on r. 
Thus Eg. <BD|1 leads to 

p=I[0 2 + 2 y(0)], p =I[0 2 -2K(0)]. (22) 

In terms of the scalar field, the integrated total energy density of the shell 
at some time is 

m = 2nR 2 [<fi 2 + 2V ((/>)]. (23) 
Using Eqs.([22D and ([23]) in Eq. flH) . we get 

This is the Klien-Gordon (KG) equation, D0 + ^ = 0, in coordinate system 
of the shell metric (4). In terms of the scalar field, the effective potential is 



Veff(R) = 1- 



M+-M_ \ 2 fQY (M + + M_) 



2vr J R 2 (0 2 + 2^/(0))/ \RJ R 
- ^R(0 2 + 2F(0))] 2 . (25) 

Now we solve the KG equation ([24]) and equation of motion f[T71) (with 
Eq.( [25i) ) simultaneously for 0(r) and R(r). In this case, the exact solution is 
not possible. We solve these equations numerically by assuming the following 
initial conditions: 0(1) = 0.19 and 0(1) = R(l) = 1. The graphs of these 
equations for the set of initial data are shown in Figure 2. The left graph 
shows the behavior of the shell radius R in which upper and lower curves 
represent the expanding and collapsing shell, respectively. The right graph 
is the behavior of scalar field whose upper and lower curves represent the 
collapsing and expanding shell, respectively. In case of collapse (upper curve), 
scalar field density goes on increasing while in case of expansion (lower 
curve), this comes to a point on r-axis implying that scalar field decays to 
zero value in this case. 

3 Massless Scalar Field 



A scalar field becomes massless, when scalar potential, V(<p), is zero. In this 

^0 = whose solution is = j^, 



case, the KG equation reduces to + ^0 = whose solution is w - - L> 
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Figure 3: The left graph shows the shell radius for massless scalar field 
case (Eq.( l26i) ). The right graph is the effective potential for massless scalar 
field (Eq.(l28l)) with Q = 1, keeping all the remaining parameters and initial 
conditions fixed as in Figures 1,2. 




Figure 4: This figure describes the behavior of effective potential (Eq.( )28p ). 
Both graphs correspond to varying M + and M_, keeping the remaining pa- 
rameters fixed as in previous cases. 



where fl is an integration constant. Thus the equation of motion (fTTj) with 
Eq.( |25l) takes the form 

M+-M\ 2 d4 (Q\ 2 (M+ + M_) -rr 2 n 4 



We define the following two parameters: 

\M\=M + -M„, M= M * + M -. 



Using these in Eq.( l26|) . it follows that 

R 2 + V eff = 0, (27) 
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Figure 5: The left graph describes the effective potential for massless scalar 
field with different values of Q, keeping all the remaining parameters fixed 
as in previous cases. The right graph represents the behavior of the massless 
scalar field shell for different values of the charge Q. 



where 

K// - 1_ U^J R + {r) --R--W (28) 

For the initial data of the shell, the left graph in Figure 3 shows the 
increase and decrease in shell radius implying the expansion and collapse of 
the massless scalar field shell. Thus a massless scalar shell may expand or 
collapse depending on the sign of velocity (i.e., R) of the shell with respect 
to stationary observer. The behavior of the potential depends on the number 
of roots of the potential. If there is no root then the scalar field shell either 
expands indefinitely or collapses to a zero size from some finite value. If there 
is one non- degenerate root then the shell expands to infinity or contracts to 
some finite size. For one degenerate root, the shell will be in an unstable 
equilibrium or collapses to form a black hole or naked singularity |30j . 

The graphical representation of the effective potential with fixed param- 
eteric values of the model is shown in Figures 4, 5. Both graphs for varying 
M + and M_ in Figure 4 and the left graph in Figure 5 show that the effective 
potential diverges for initial values of R and then V e ff — > — oo as R — » oo. In 
these cases, the shell expands to infinity or collapses to zero size. The right 
graph in Figure 5 shows that the effective potential has one root and there 
occurs unstable situation, after which potential diverges negatively and shell 
expands or collapses. The cases in which collapse occurs, the shell collapses 
to zero size by forming a curvature singularity at which intrinsic Ricci scalar 
of the shell, = -j^(2RR + R 2 ), diverges. 
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Figure 6: The behavior of effective potential for massive scalar field is shown 
for fixed parameters as well with varying charge parameter. 

4 Massive Scalar Field 



In this case, we discuss the motion of a scalar field for which potential term, 
V(<j)) is determined by taking p as an explicit function of R. From Eq. (122j) . 

we get 

12 „ l ( „ ^ (2Q) 



p + p, V ((/>) = ~(p- p). 



Also, from Eqs. ffT2]) and ( T3~3j) . we get 



dp 2 

dR + R iP + p) 



0. 



(30) 



Here we use p as an explicit function of R, [29] i.e., p = p e kR , where p 
and k are constants. Inserting this value of p in Eq. fl30|) . it follows that 



P= R- 



X , 2(1 + kR)p e 



-kR 



k 2 R 2 



(31) 



where x is constant of integration. Notice that the above equations satisfy 
the conservation equation (TH|) . Further, applying the values of p and p in 
Eq.fl2H]), we get 



v{4>) 



x 

2R 2 



-kR 



2 

-kR 



1 



1 + 



2(1 + kR) 
k 2 R? 
2(1 + kR) 
k 2 R 2 



(32) 
(33) 



11 



These equations satisfy the KG equation Q2H). Using Eqs.fl22J-(|32J in (J2SD, 
we have 



( 



M. 



M. 




2 



2 



(M+ + M_) 
R 



^e//(i2) = 1 



+ 



m 



where 



87rp e 



-kR 



m = AnR 2 p = Anx + 



(1 + kR). 



(35) 



The behavior of effective potential for massive scalar field shell is shown 
in Figure 6. The left graph is effective potential for massive scalar field 
(Eq.(l34j)) for k — 1, x = 3,£>o = 1 an d remaining parameters are fixed as in 
the massless scalar field case. This implies that V e ff — > — oo as R — > 0, the 
massive shell collapses to zero size forming a curvature singularity. The right 
graph represents effective potential for massive scalar field shell for different 
values of Q. There appear oscillations in the system. There exist such values 
of charge parameter for which scalar field shell executes an oscillatory motion. 
The oscillations occur at two points where V e ff cuts the horizontal axis at 
more than one point. The values of R for which V e ff = are shown in right 
graph of Figure 6 yielding zero velocity. This implies that the shell stops 
for a moment and then expands or collapses. During the collapsing phase at 
minimum values of the radius, the tangential pressure reaches its maximum 
values while during the expansion, minimum pressure occurs at maximum 
radius. In this way, scalar field shell performs the oscillatory motion. The 
values of R for which V e ff = 0, and intrinsic curvature of the shell is finite, 
are bouncing points after bounce the shell either expands or collapse. 

5 Discussion 

In this paper, we have examined the dynamical behavior of the scalar field 
thin shell. Using the Israel formalism, the equations of motion have been for- 
mulated by taking the internal and external regions to the boundary surface 
as RN solution. The equations of motion are originally derived for perfect 
fluid and then are written in terms of scalar field. The complete dynamics of 
the thin shell is described by the equation of motion (JT7J) and the KG equa- 
tion ( 124")) . The exact solution of these equations cannot be found, however 
can be solved numerically. Firstly, we have solved these equations by using 
the scalar field potential as quadratic potential. This solution is shown in 
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Figures 1, 2 which represents the collapsing and expanding scalar field shell. 
It has been found that scalar field decays out in the case of expansion while 
it grows in the case of collapse. 

To analyze further, we have taken massless scalar field and a massive 
scalar field with p as an explicit function of R. In the case of massless scalar 
field, it has been found that the shell radius is an increasing or decreasing 
function of the proper time implying that the shell either collapses or ex- 
pands. In this case, the effective potential (Figures 3-5) predicts that shell 
can collapse to zero size by forming a curvature singularity or can expand to 
infinity. For the massive scalar field case with p as an explicit function of R, 
we have evaluated the scalar field potential instead of taking it as quadratic. 
It has been found that the shell radius behaves like the massless scalar field 
and the effective potential (Figure 6) diverges negatively as the radius of the 
shell approaches to zero and oscillating behavior is noted in this case. This 
indicates that the shell collapses to zero size forming a curvature singularity. 
There is also bouncing behavior of the shell in this case. In the both (mass- 
less and massive scalar field) cases, when shell collapses, the edge of the shell 
coincides with the horizons of the interior black hole. 

The results can be summarized as follows. We have found that there 
are three possible phases (expanding, collapsing and oscillating (bouncing)) 
during the dynamics of the scalar field in the present configuration. When 
shell expands, it continues expanding forever with constant velocity as the 
boundary surface is described by the spatially homogenous spacetime. In 
case of collapse, a shell collapses to zero size forming a curvature singularity 
at R = 0, where intrinsic Ricci scalar R^^ = — -j^(2RR + R 2 ), diverges. Also, 
the turning (bouncing) points occur when V e ff(R) = at more than one value 
of R. In this case, the oscillations occur between two points where V e ff = 0. 
The values of R for which V e ff — (right graph in Figure 6), one gets zero 
velocity. This implies that the shell stops suddenly and then expands or 
collapses. During the collapse at R = R m in, the tangential pressure reaches 
to its maximal values while during expansion minimal pressure occurs R = 
R m ax- In this way, the scalar field shell performs the oscillatory motion. 

It would be interesting to extend this work for more generic geometry 
or using the polytropic equation of state to check the validity of cosmic 
censorship hypothesis. 
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